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ABSTRACT 

One  dimensional  Rectangular,  Barlett,  Cosine,  Hann, 
Hamming,  Blackman  and  Kaiser  window  are  extended  to  two 
dimensions  using  circular  symmetry  and  the  Hankie  transform 
in  integral  form.   It  is  shown  that  this  transform  cannot 
generally  be  expressed  in  closed  form.   Comparisons  of  the 
above  two  dimensional  windows  are  then  made  on  the  basis 
of  two  dimensional  extensions  of  the  criteria  of  equivalent 
noise  bandwidth,  coherent  gain,  and  the  so  called  trade-off 
parameters . 
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I.   INTRODUCTION 

Windowing  is  essential  to  signal  processing  and  to  digi- 
tal filter  design.   For  example,  the  weighting  functions  of 
the  window  are  used  in  harmonic  analysis  of  signals  [1]  and 
in  design  techniques  for  linear-phase  digital  FIR  filter  [2]. 

In  fact,  if  an  infinite  sequence  of  numbers  is  truncated, 
then  its  Fourier  transform  will  be  subject  to  the  well-known 
Gibb's  phenomenon;  which  will  cause  over  shoot  and  ripples 
before  and  after  any  point  of  discontinuity  in  the  frequency 
domain.   The  overshoot  and  ripples  will  not  disappear  no 
matter  how  long  the  sequence  is . 

To  overcome  this  undesirable  effect  in  the  frequency  do- 
main, the  finite  sequence  should  be  modified  using  a  weighting 
function,  called  a  window.   For  example,  if  w(n)  is  a  weight- 
ing sequence  and  h(n)  is  an  infinite  sequence  of  number,  then 
the  windowed  finite  sequence  is: 

h  (n)  =  w(n)  .  h(n)  (1.1) 

w 

where : 

f(n)         |n|<N 


w(n)  =    {   0  elsewhere 


(1.2) 


In  this  example  the  weighting  functions  performs  both 
truncation  and  modification  of  the  signal  at  the  same  time. 
Multiplication  of  two  sequences  in  time  domain  is  equivalent 


to  the  convolution  of  their  transforms  in  the  frequency 
domain.   Consequently  the  windowing  method  has  the  effect 
of  smoothing  out  the  rough  points  -  moderating  the  overshoot 
and  ripples  -  in  the  original  frequency  response. 

The  technique  of  windowing  is  shown  graphically  in 
Fig.  1.1. 

There  are  many  window  functions  available  in  many  pieces 
of  literature.  [1,  3]   According  to  the  way  in  which  the 
window  function  is  formed,  they  can  be  separated  into  two 
categories;  basic  windows  and  constructed  windows.   A  con- 
structed window  is  formed  from  a  product,  a  sum,  or  a  con- 
volution of  two  or  more  basic  windows. 

Seven  best-known  basic  windows  are  listed  in  Table  1.1, 
and  they  are:   Rectangular,  Barlett,  Cosine,  Hann,  Hamming, 
Blackman  and  Kaiser  window.   Most  of  the  windows  are  named 
after  individuals.   These  seven  basic  windows  will  be  dis- 
cussed thoroughly  in  later  chapters. 

Although  Huang  [4]  shows  that  a  good  two  dimensional 
window  can  be  derived  from  one  dimensional  window,  most 
windows  found  in  the  literature  are  given  in  only  one-di- 
mensional form.   It  is  the  purpose  of  this  thesis  to  present 
and  discuss  two-dimensional  windows  derived  from  those  one 
dimensional  form,  using  Huang's  theorem. 

The  easiest  way  to  compare  various  windows  is  to  create 
some  appropriate  figures  of  merit  for  windows.   This  thesis 
presents  those  figures  of  merit  that  are  used  to  compare  the 
given  window  functions. 
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A.   TECHNIQUE  OF  WINDOWING  [2] 

A  graphical  interpretation  of  windowing  is  given  in 
Fig.  1.1  (a  non-recursive  filter  design  using  the  window 
method) . 

The  first  figure  shows  an  infinite  impulse  response  of 
a  digital  filter,  h(n)  ,  and  its  periodic  frequency  response 
H(w)  . 

The  next  figure  shows  a  rectangular  window  r(n) ,  with 
its  frequency  response,  R(w) . 

The  third  figure  consists  of  multiplication  of  h(n)  and 
r(n) ,  which  produces  a  truncation  of  the  filter  impulse 
response,  h  (n)  .   Multiplication  in  the  time  domain  is  equi- 
valent to  convolution  in  the  frequency  domain.   The  convolu- 
tion of  H(w)  and  R(w)  is  H  (w)  containing  overshoot  and 
ripples  in  its  frequency  response  of  the  corresponding  digi- 
tal filter  as  shown  in  the  figure. 

Truncation  of  impulse  response  is  required  for  a  non- 
recursive  digital  filter  design,  since  the  filter  produces 
an  output  which  depends  only  upon  the  present  and  a  finite 
number  of  past  inputs . 

The  fourth  figure  shows  the  non-rectangular  window  w(n) , 

and  its  periodic  Fourier  Transform  W(w) .   The  purpose  of 

this  window  is  to  alter  the  truncated  filters  impulse  response, 

h. (n),  slightly  in  order  to  minimize  the  overshoot  and  ripples 

in  its  frequency  response  H  (w) . 

The  last  figure  shows  the  finite  impulse  response  h  (n) , 

w 

which  consists  of  multiplication  of  h  (n)  and  w(n) .   The 
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frequency  response  of  h  (n)  is  produced  by  convolving  H  (w) 

with  W(w) ,  resulting  in  a  smooth  frequency  response  H  (w) . 

w 

The  technique  of  windowing  applied  to  a  filter  design, 
results  in  an  improved  filter  with  a  smooth  frequency  response. 

B.   THESIS  ORGANIZATION 

The  following  discussion  presents  the  organization  of 
this  thesis. 

As  mentioned  before,  the  best  way  to  discuss  windows  is 
to  define  parameters  which  can  be  used  to  characterize  a 
window.   Logically,  the  window  figures  of  merit  are  first 
discussed  in  the  next  chapter,  Chapter  II. 

Chapter  III  gives  an  overview  of  one-dimensional  windows 
and  specifies  their  characteristics,  using  parameters  defined 
in  Chapter  II.   Table  III.l,  given  at  the  end  of  this  chapter, 
tabulates  one-dimensional  windows  and  their  parameters. 

The  transformation  from  one-dimensional  to  two-dimensional 
is  given  in  Chapter  IV,  where  we  present  three  major  discus- 
sions.  First  we  present  the  expansion  of  the  Fourier  Trans- 
form from  one-dimension  to  two-dimensions.   Then  we  present 
how  Huang  [4]  expanded  a  good  one-dimensional  window  into  a 
good  two-dimensional  window.   Finally,  the  definition  of  one- 
dimensional  figures  of  merit  are  expanded  into  two-dimensional 
window  figures  of  merit. 

Chapter  V  presents  two-dimensional  windows  in  detail,  and 
specifies  their  parameters  as  defined  in  Chapter  IV.   Table 
V.l,  given  at  the  end  of  this  chapter,  tabulates  two-dimen- 
sional windows  and  their  parameters. 
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Figure    1.1 

Application  of  Window  Function  To 
FIR  Filter  Design 
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Two-dimensional  Fourier  Transform  of  window  functions 
are  hard  to  present  in  a  closed  form.   Fortunately  a  closed 
form  Fourier  Transform  is  not  required  in  studying  windows, 
so  two-dimensional  Fourier  Transform  of  the  windows  are 
presented  in  integral  form.   The  effort  to  present  a  two- 
dimensional  Fourier  Transform  of  a  window  in  a  closed  form 
is  given  in  the  Appendix. 


li+ 


TABLE    1.1 
Seven   Basic   One-Dimensional  Window   Functions 


No.        Window  Function 


1 

Rectangu 

1 

ar 

W(x)  =  1.0 

2 

Barlett 

W(x)  =  1.0 

+  M 

a 

3 

Cosine 

W(x)  =  cos 

TT 

IE   X 

I x  I  _la 
I x  l_5a 

h/      \  11  TT  11^ 

W(x)    =   75-  +   7T  cos   —  x  x    <a 

2         2  a  11— 

5  Hamming  W(x)    =    0.54    +    0.46    cos   —  x         |x|<a 

ci 

6  Blackman  W(x)    =    0.42+0.5    cos   ^  x    +    0.08    cos    -  x 

2  a 


x  I  <a 


I 


{3    /    1.0-(^)2} 


Kaiser  W(x)    =   — = — t-^-t lxlj;a 
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II.   WINDOW  FIGURES  OF  MERIT  [1,  3] 

Since  windows  are  used  in  frequency  response  analysis  to 
reduce  undesirable  effects  -  overshoots  and  ripples  in  the 
frequency  domain  -  it  is  very  important  to  find  the  para- 
meters or  figures  of  merit  of  the  windows  that  can  be  used 
to  identify  their  characteristics  and  to  indicate  their  per- 
formance . 

There  are  many  convenient  measures  that  can  be  used  as 
parameters . 

The  common  desirable  characteristics  of  a  window  are: 

1.  A  narrow  main  lobe  width  of  the  frequency  response 
of  the  window,  containing  as  much  of  the  total  energy  as 
possible . 

2.  Side  lobes  of  the  frequency  response  decreasing  in 
energy  rapidly  as  frequency  increases. 

These  two  common  characteristics  certainly  can  be  used  as 
parameters  to  identify  the  window's  performance. 

Harris  [1]  and  Geckinli  [3]  define  one-dimensional  window 
performance  parameters  which  will  be  repeated  here.   These 
parameters  are  used  throughout  this  thesis  and  will  be  expan- 
ded for  the  two-dimensional  case  in  later  chapters.   The 
parameters  are:   equivalent  noise  BW ,  coherent  gain,  and 
trade-off  parameters. 
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A.   EQUIVALENT  NOISE  BANDWIDTH  (ENBW) 

The  ENBW  parameter  is  the  width  of  a  rectangular  filter 
in  the  frequency  domain  with  the  same  amplitude  as  the  peak 
power  gain  of  the  window,  and  with  a  width  such  that  the 
total  noise  power  in  the  rectangular  filter  is  equal  to  the 
total  window  noise  power. 

A  graphical  explanation  of  the  parameter  is  shown  in 
Fig.  2.1.   The  first  figure  shows  a  window  frequency  response 
which  consists  of  a  main  lobe  and  side  lobes.   The  second 
figure  shows  a  discrete  frequency  component  of  a  signal, 
together  with  white  noise.   The  third  figure  shows  the  result 
of  windowing. 

The  window's  frequency  response  convolves  with  the  dis- 
crete component  of  the  signal  plus  noise.   In  this  case  the 
window  acts  as  a  filter,  gathering  information  over  its 
bandwidth.   The  last  figure  represents  the  area  of  the  windowed 
discrete  frequency  response  of  the  signal  plus  noise,  at  dis- 
crete frequency  W  . 
^       o 

It  is  clear  from  the  figure  that  to  minimize  the  accumu- 
lated noise  collected  by  the  window,  it  is  required  that  the 
window  has  a  narrow  bandwidth. 

A  convenient  measure  of  window's  bandwidth  is  called 
equivalent  noise  bandwidth. 

Noise  power  accumulated  by  the  window  is  defined  as: 

N 
Noise  Power  =  -#  /    |W(w)|    dw  (2.1) 

where  N   =  noise  power  per  unit  BW. 
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Figure  2.1 

Graphical  Interpretation  Of  Convolution  Of  A  Window 
And  Spectral  Amplitude  Of  A  Signal,  At  w=wq 
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Parseval's  theorem  allows  equation  (2.1)  to  be  computed  as: 

N     T 
Noise  Power  =  -£        f      WZ(t)  dt  (2,2) 

The  integral  in  (2.1)  and  (2.2)  represents  the  area  under 

the  window  curve.   If  the  area  is  divided  by  height,  then 

the  result  is  the  width  of  a  rectangle  that  has  the  same  area 

with  the  window. 

The  peak  gain  of  signal  occurs  at  w=0,  and  is  defined  as: 

T 
Peak  signal  gain  =  w(0)  =   /  w(t)  dt  (2.3) 

-T 

2  2 

Peak  power  gain   =  w  (0)  =  [  /  w(t)  dt]        (2.4) 

-T 

Then  the  ENBW  is  equal  to  the  noise-power  divided  by  peak 

power: 

T 

/  w2(t)  dt 

ENBW   =   ^-^ 2  (2-5) 

[   /  w(t)  dt] 

-T 

The  above  expression  is  a  normalized  ENBW,  after  it  is 

.  .  N 

divided  by  factor  _o  .   This  parameter  will  be  used  as  a 

2T 
window  figure  of  merit. 

B.   COHERENT  GAIN 

Coherent  gain  is  best  defined  using  a  sinusoidal  signal. 
In  Fig.  2.1  a  component  of  signal  frequency  response  at  a 
discrete  frequency  is  shown.   Convolving  this  signal  with  a 
window's  frequency  response  can  be  thought  as  passing  the 
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W(w) 


Peak  Power  Gain  =  |W(o)| 


Equivalent  Noise  BW 


w 


Figure  2.2 
Equivalent  Noise  BW  Of  A  Window 
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signal  through  a  filter.   Using  this  view,  the  coherent  gain 
of  a  window  can  be  examined  by  passing  a  sinusoidal  signal 
through  a  filter. 

Let  the  sinusoidal  signal  be  defined  as: 

f(x)  =  Aejwx   +   n(x)  (2.6) 

where  n(x)  =  white  noise. 

The  filter  output  on  the  windowed  signal  is  given  as: 

f  (x)  =  A  w(x)e^WX   +   n(x)  w(x)   -   a<x<a     (2.7) 
w 

The  frequency  response  of  the  windowed  signal  only  (less  noise) 

is  given  as: 

a 
F  (f)   =   /A  W(x)e^wx  e"]WX  dx   = 
w        -a 


a 
A   /   W(x)  dx  (2.8) 

-a 


The  windowed  frequency  response  is  proportional  to  the  input 
amplitude  A.   The  proportionality  factor  is  the  integral  of 
the  window  function. 

For  a  rectangular  window  this  factor  is  2a.   For  any  other 
window  the  factor  is  reduced  due  to  the  window  smoothly  going 
to  zero  near  the  boundaries. 

The  coherent  gain  is  defined  as  the  proportinality  factor 
divided  by  factor  2a,  the  extension  of  the  window  (normalized). 

1    a 
COH.GAIN   =   -g=-    /   W(x)  dx  (2.9) 

c  a 

-a 
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The  coherent  gain  is  the  second  parameter  that  will  be  used 
as  window  figure  of  merit . 

C.   TRADE-OFF  PARAMETERS 

Trade-off  parameters  are  best  explained  using  Fig.  2.3. 
The  top  figure  shows  the  Hamming  window  in  the  time  or  spa- 
tial domain,  and  the  bottom  figure  shows  its  Fourier  trans- 
form, which  is  given  in  log  frequency  versus  normalized  log 
magnitude  (20  log  W(f)/W(0)). 

The  window's  characteristics  in  frequency  domain  are 
largely  determined  using  these  trade-off  parameters. 

From  Fig.  2.3(b)  these  parameters  are: 

1.  a    :•   highest  side  lobe  level  in  dB. 

2.  a.         :   side  lobe  level  at  a  higher  frequency, 

where  asymptotic  assumption  is  valid  for 
the  window  behavior. 
3  .   b    :   frequency  at  which  main  lobe  drops  to  the 

peak  side  lobe. 
M- .   d    :   side  lobe  fall-off  rate  in  dB/octave. 
Clearly  from  the  above  definitions ,  the  smaller  the  value 
for  b,  the  better  the  window  performance;   the  smaller  the 
value  for  a-, ,  the  smaller  the  leakage  through  the  near  side 
lobes ;  the  smaller  the  value  for  a? ,  and  the  larger  the  value 
of  d,  the  smaller  the  leadage  through  the  far  side  lobes. 
These  parameters  are  expanded  for  two-dimensional  case  in 
a  later  chapter. 
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Figure  2.3 

Plot  Of  Hamming  Window,  Is  Used  To  Illustrate  The 
Definition  Of  Trade-Off  Parameters  a, ,  a? ,  b  and  d. 
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III.   ONE  DIMENTIONAL  WINDOW  FUNCTIONS 
[1,  3,  5,  6,  7] 

The  seven  common  window  functions  listed  in  Table  1.1, 
are  reviewed  in  detail  in  this  chapter.   Figures  of  merit 
for  each  window  are  presented  using  the  definitions  given  in 
the  preceding  chapter.   The  window  functions  treated  here  are 
in  continuous  form  rather  than  in  discrete  form. 

Without  loss  of  generality,  let  the  window  function, 
w(x) ,  be  unity  at  the  origin  and  limited  to  the  interval 
|x|_<  a.   Presented  in  mathematical  form  this  means: 

00 

w(o)   =  .  /   W(f)  df  =  1  (3.1) 

—  00 

w(x)   =   0  for  |x|>a  (3.2) 

The  Fourier  Transform  of  the  above  window  is  given  as : 

a 
W(f)   =   /   w(x)  e"]wx  dx     |x|_<a  (3.3) 

-a 

The  transforms  of  the  windows  generally  consists  of  a  main 
lobe  with  a  large  peak  amplitude  and  side  lobe  with  compara- 
tively small  peak  amplitudes.   To  better  display  the  presence 
of  the  side  lobes,  logarithmis  scales  are  used  for  plotting 
in  place  of  linear  scales . 
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A.   RECTANGULAR  WINDOW  [1] 

The  rectangular  window  is  commonly  called  a  boxcar  or 
Dirichlet  window.   Its  main  purpose  is  to  truncate  the 
length  of  a  function  to  a  finite  extent,  without  affecting 
its  amplitude. 

The  expression  for  this  window  is  given  as : 

w(x)   =   1.0  l.xlla  (3.4) 

The  Fourier  Transform  of  the  window  is  given  as : 


11/ *\  1      r  -1WX   , 

W(f)   =  tt-        I      e  J    dx 

l.  a 

-a 


sinaw  f       »  / o  r  % 

=   =   sine  (aw)  (3.5) 

aw 

The  transform  is  seen  to  be  a  Dirichlet  kernel.   In  this 
report  the  extent  of  the  window  is  taken  as  a=0.5,  and  the 
ENBW  is  given  as : 


0.5 
/    dx 

ENBW   =   ^-f o    =   1-0 

[   /   dx  ] 

-0.5 


The  coherent  gain  is  given  as: 

0.5 
COH.GAIN  =   /    dx   =   1.0 
-0.5 
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Plots  of  the  window  function  and  its  tranform  are  given 
in  Fig.  3.1(a)  and  (b)  respectively.  The  main  lobe  width 
(the  width  between  origin  and  the  first  zero  crossing)  is 
1.0,  and  the  first  side  lobe  peak  is  approximately  13  dB 
down  from  the  main  lobe  peak.   The  side  lobe  fall  off  is 
6.0  dB/octave. 

B.   BARTLETT  WINDOW 

The  Bartlett  window  is  also  called  Fejer  or  triangular 
window,  and  can  be  thought  of  as  a  convolution  of  rectangu- 
lar window  with  itself.   This  window  is  the  simpliest  of  all 
windows  that  changes  the  data  values  by  multiplication  by  a 
number  other  than  unity.   The  mathematical  expression  of  this 
window  is  given  as: 

w(x)  =  1.0  -  -L^-L  |x|  <a  (3.6) 

a  — _ 

The  transform  of  this  window  is  given  as: 

-.a         II' 

W(f)  =   -     f    (1.0  -  -L2JL)e  ]WX  dx   = 

a  a 

-a 

..   o  .         t   a 

1    r        /-in     Xn-IWX         1    r  , ,   n     Xs-jWX, 

-    —      I      (1.0  +  — )e  J    dx  +  —   /  (1.0  -  — )e  J    dx   - 
a  a  a  a 

-a  o 

,o  t   o       .        ,  a 

1   ,   -iwx       1    f  -iwx  ,     1  ,   -dwx  , 

=  -/eJ  dx  +  ^5-/xeJ  dx  +  -  /  e  J    dx  - 
a                ^                  a 

a  a   -a  o 

-i    a 

1    ,    -iwx  , 
-  — j     J    x  e  J    dx   = 

2a   o 

i   p-jwx  o         o    dp-3wx    1   e-3wx  a 
a    -3w  J       2         -3w     a    -Dw   ^ 
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20    login|W(f )/W(0) | 

1U  b    =    0.81 


log    f 


(b) 


Figure  3.1 
(a)  Rectangular  Window,  (b)  Spectral  Window 
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.  a  ,  -jwx 

1  ,  de  J 

-  — o  /  x  = 

2         -iw 

a  o       J 


1   -|_eJaw     l   ve"^wx  ?     °   e"3WXdx,    1  e-:iaw-l 
a    -jw      72     -jw   '  -]w      a   -jw 

d  —  d.       — a. 


,      -iwx  a  a  -iwx 

—75-  { » -  /   = dx} 

2     -iw  -iw 

a       J    o  o  J 


,   law   -naw     ,    law   -law  ,  o  -iwx 

l,eJ   -e  J   v    1  ,eJ   -e  J   »  1  ,  e  J    , 

— ( )  -  —  ( . )  -  — ■*■  j  ■. dx  + 

a               a       iw  2  -iw 

J  a  -a    J 


,  a  - i wx 

+  — yr  J  . dx 

2  -iw 

a  o    J 


.    .   2  o  ,  .  2    .    a 

.  1  (_i_)  e-3wx  !  +  1  (_1__  e-]wx  ! 

2  -iw  '  2  -iw 

a    J  a  a  J  o 


+   ^_   (1    .    e  +  i^   .    e-3aw   +    1}    =   _zl_   (e:aw   _    2    +    e-jaw} 
aw  aw 

2    jaw/ 2   -jaw/ 2  2 

^   (- =53 >      • 


2 
,  sin  aw/ 2  >.     ,     .  ,^2  ,.  -, 

( 7-0 )   =  (sine  aw/2)  (3.7) 

aw/  z 


The  transform  is  seen  to  be  the  squared  of  Dirichlet  Kernel. 
If  the  extent  of  the  window  is  taken  as  a=0. 5,  then  the 
parameters  ENBW  and  COH.GAIN  are: 
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0.5  2 

/   (1-  q^)   dx 

ENBW  =  ^s : 2    =  krt  =  x-33 

[  /  (i  -  JL  dx  ] 

-0.5       °'5 


0.5 
COH.GAIN  =   /    (1  -  —pO  dx  =  °-5 
-0.5        U,i) 


The  window  function  and  its  transform  are  shown  in  Fig. 
3.2(a)  and  (b)  respectively.   Its  main  lobe  width  is  2.0, 
twice  that  of  the  rectangular  window,  and  the  side  lobe  fall- 
off  is  12  dB/octave. 

C.   COSINE  WINDOW 

Cosine  window  is  actually  one  member  of  a  family  window 
with  a  general  form  given  as: 

w(x)  =  cos   (75— x)  Ix.l  <a  (3.8) 

2a  11  — 

Type  of  this  window  depends  upon  parameter  a,  which  usually 
has  a  value  between  1  and  4.   Of  particular  interest  in  this 
family  window  is  Hann  window,  where  a=2,  which  will  be  con- 
sidered later. 

The  cosine  window  is  defined  with  cc  =  l,  as: 

w(x)   =   cos  (-*— x)  lxl<a-  (3.9) 

z  a  — 

The  Fourier  Transform  of  this  window  is  given  as: 


i   a 

w(f)  =  -7T—     J      cos  -7T- x  e  J    dx 
z  a         zd 
-a 
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20    loglQ|W(f)/WCO)|  g 


llifcl 


(a)    Barle 


(b) 

Figure    3 . 2 
tt   Window,    (b)    Spectral   Window 


1   /   e-3(w-^)x  dx  ♦  i-  /   e"3(w+5S)x  dx 
Za  za 

-a  -a 


7T  TT 

,  sin  a(w— s— )     -■  sin  a(w+7r— ) 
1  2a    ,  1  2a 


a(w-7i°  a(2+2¥) 


■s-  sine  a(w-~— )  +  -=-  sine  a(w+-~— )  (3.10) 

Z  Z  3.  Z  ZcL 


The  transform  of  this  window  is  a  summation  of  two  shifted 
Dirichlet  Kernel. 

If  the  extend  of  the  window  is  taken  as  a=0.5,  then  the 
parameters  ENBW  and  COH.GAIN  are  given  as: 


0.5  2 

/         (costtx)      dx 

enbw  .   ^5 2    ■   rw  =    x-23 

[      /        costtx   dx] 
-0.5 


0.5  9 

COH.GAIN      =      /      costtx   dx      =      -      =       0.64 


-0.5 


TT 


Plots  of  the  window  and  its  transform  are  given  in  Fig. 
3.3(a)  and  (b)  respectively.   From  the  transform  plot,  the 
main  lobe  width  is  1.5  times  that  of  the  rectangular  window, 
The  first  side  lobe  peak  is  approximately  2  3  dB  down  from 
the  main  lobe  peak,  while  the  side  lobe  fall  off  rate  is 
12  dB/octave. 
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Figure  3  .  3 
(a)  Cosine  Window,  (b)  Spectral  Window 
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D.   HANN  WINDOW 

The  Hann  window  is  one  member  of  the  family  windows  given 

ex  •  ... 

by  cosx,  where  a=2.   The  window  equation  is  written  as: 


/     \  2ttx  11  TTX  t     i  ^  ,on^N 

w(x)       =      cos    o—      =      o    +    o    cos   —  x    <a  (3.11) 

2a  2         2  a  '     '— 


The  Fourier  Transform  of  this  window  is  given  as 


WCf )      =      ~     /      (h   +   h   cos-  x)    e-;]WX   dx 
za  a 

— a 


,        a  ,  a 

1            -  -]WX                        1  r                      TT            -JWX 

r—      f  e    J         dx    +   t—  /  cos   — x   e    J         dx 

4  a  4a                       a 

-a  -a 


TT  TT 

-,       .  ,    sin   a(w — )  ,    sin   a(2+— ) 

lsin_aw      +1  a_  +1 a_  (  } 

2       -    aw  4  ,        TTx  4  f      .TTv 

a(w — )  a(w+— ) 

a  a 


The  Transform  of  this  window  consists  of  a  summation  of 
three  Dirichlet  Kernels ,  one  at  the  origin  and  the  two  others 

TT 

shifted  —  radians  from  the  origin, 
a  & 

The  three  terms  on  the  transform  of  this  window  are 
shown  in  Fig.  3.4,  where  the  shifted  Dirichlet  Kernels  are 
located  at  the  first  zeroes  of  the  center  kernel.   The  peak 
amplitudes  of  the  shifted  kernels  are  half  that  of  the  center 
kernel;  also  the  side  lobes  of  the  translated  kernels  are 
about  half  the  size  of,  and  of  opposite  phase  from  those  of 
the  center  kernel.   The  summation  of  the  three  side  lobes  in 
phase  opposition,  tends  to  cancel  the  overall  side  lobe  structure 
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i  W(f) 


0.25  sine  a(w+— ) 


0 . 5  sine  aw 


i  \     0.25  sine  a(w--) 


Figure  3.4 

Transform  Of  Hann  Window  As  A  Sum 
Of  Three  Dirichlet  Kernel 
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The  cancellation  of  the  side  lobes  structure  suggests 
a  constructive  technique  to  define  new  windows .   The  best- 
known  windows  with  this  type  of  side  lobes  cancellation  are 
Hamming  and  Blackman  windows ,  which  are  given  in  the  next 
section . 

If  the  extent  of  the  Hann  window  is  taken  as  a=0.5,  the 
parameters  ENBW  and  COH.GAIN  are  given  by: 


0.5  2 

/  (h+h    cosp— f-)   dx 
-OS  0375 

ENBW  =  -Hrs : =  to  =  1-5 

[    /  (h+h    cos^x)  dx]Z 
-0.5  U*5 


0.5 
COH.GAIN  =     f         (h+h    cos  J^r    )  dx  =  0.5 
-0.5  U,D 


The  window  and  its  Fourier  Transform  are  plotted  in  Fig 
3.5.   The  main  lobe  width  is  2.0,  twice  that  of  the  rectan- 
gular window  and  the  first  side  lobe  level  is  23  dB  down 
from  the  main  lobe.   The  side  lobe  level  fall  off  rate  is 
18  dB/octave. 
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Figure  3 . 5 
(a)  Hann  Window,  (b)  Spectral  Window 
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E.   HAMMING  WINDOW  [5] 

Hamming  window  can  be  thought  an  an  improvement  of  Hann 
window,  in  a  sense  that  the  first  side  lobe  cancellation 
is  even  greater.   From  Fig.  3.M-,  the  summation  of  the  three 
side  lobes  produce  an  inexact  cancellation.   The  cancellation 
can  be  made  more  pronounced  by  adjusting  relative  size  of 
each  term. 

Exact  cancellation  of  the  side  lobes  was  first  shown  by 
Hamming  using  a  window  given  in  general  form  as: 


IT  i   i 

w(x)  =  a  +  (1-a)  cos  —  x  | x | <a  (3.13) 

a  ""■ 


where  a  is  a  coefficient  that  will  be  defined  later.   The 
Fourier  transform  of  the.  window  is  given  as: 


7T 

W(f)  =  asinc  aw  +  0.5(l-oc)sinc  a(w — ) 

a 

+  0.5(l-a)sinc  a(w+J)  (3.14) 

a 


From  (3.14)  the  best  cancellation  of  the  first  side  lobe 

25  25 

occurs  when  a  =  ^r  ,  {-^-~    =  0.543478261).   If  a  is  selected  as 

26  26 

25 

0.54,  an  approximation  of  y^-,  then  an  improvement  in  cancel- 
lation of  the  side  lobes  occurs.   The  value  a=0.54,  gives 
rise  to  a  window  called  the  Hamming  window,  which  is  written 
as : 


w(x)  =  0.54  +  0.46  cos  —   .         I x I  <a  (3.15) 

a  ii_ 
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The  Fourier  transform  then  is  given  as: 


W(f)  =  0.54  sine  aw  +  0.23  sine  a(w--)  + 

a 


+  0.23  sine  a(w+-)   .  (3.16) 

a 


If  the  extent  of  the  window  is  a  =  0.5,  then  the  para- 
meters ENBW  and  COH.GAIN  are  given  as: 

0.5  2 

/    (0.54+0.46cos  n-^-r-x)   dx 
-0  5  n  3974. 

ENBW  =  ^irh 7  '-  Iritis  =  i.as 

[   /   (0.54+0.46cos7^Tx)  dx] 
-0.5  U*b 


0.5  2 

/   (  0.54  +  0.  4  6cos?f^Tx)   dx 

-0  5  0  3974 

COH.GAIN  =  ^_ _  =  0^3974  =  1>36 

[   /   (O.SU  +  O.^ScosTr^-x)  dx] 
-0.5  u,b 


The  Hamming  window  and  its  transform  are  ploted  in  Fig. 
3.6(a)  and  (b)  respectively.   From  its  transform  plot,  the 
main  lobe  width  is  twice  that  of  the  rectangular  window,  and 
the  first  side  lobe  level  is  4  3  dB  down  from  the  main  lobe. 
The  side  lobe  fall  of  rate  is  6  dB/octave. 

F.   BLACKMAN  WINDOW  [6] 

Hann  and  Hamming  windows  are  examples  of  windows  con- 
structed as  the  summation  of  shifted  sinc(x)  functions. 

Blackman  suggested  a  general  rule  by  which  window 
functions  could  be  constructed  as  follows: 
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Figure  3 . 6 
(a)  Hamming  Window,  (b)  Spectral  Window 
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w(x)  =    Z   a   cos  m  ( — )  x  <a        (3.17) 

-n   m  a 

m=0 


and  subject  to  the  constraint 


£   a   =1.0  (3.18) 

m=0 


The  Fourier  transform  of  this  window  becomes: 


k       a 
W(f)  =    E  (-l)m  -#  {sine  a(w--m)  +  sine  a(w+-m)}    (3.19) 

m=0       l  a  a 


Hann  and  Hamming  windows  are  of  this  form  with  a„ ,  and  a, , 
being  non  zero.   Their  spectral  windows  are  summations  of 
three  shifted  kernels. 

Any  value  of  k  will  produce  a  new  window  with  a  summation 
of  (2k-l)  kernels  in  its  tranform,  and  k  non  zero  corficients, 
however  one  way  to  achieve  windows  with  a  narrow  main  lobe 
is  to  restrict  k  to  a  small  interger. 

Blackman  examined  this  window  with  k=3,  and  found  best 
cancellation  of  the  third  and  fourth  side  lobes  occurs  when: 


ao  "  i£m  --  °-u25  59°  71  -  °-u2 


9240 
a,  =  n  QenQ  =  0.496  560  62  =  0.50 
1     lobUo 


a2  =  wms  -  °-076  8U8  67  -  °-° 


HO 


The  window  function  which  uses  these  three  coeficients 
is  known  as  the  Blackman  window,  and  is  given  as: 


w(x)  =  0.42  +  0.5  cos—  +  0.08  cos— x     I  x  I  <a     (3.20) 


The  Fourier  transform  of  this  window  is: 


W(f)  =  0.42  sine  aw  +  0.25  sine  a(w--)  +  0.25  sine  a(w+-) 

a  a 


+  0.04  sine  a(w-— )  +  0.04  sine  a(w+— ) 

a  a 

(3.21) 

If  the  extent  of  this  window  is  a=0.5,  then  the  ENBW  and 
COH.GAIN  parameters  are  given  by: 


0.5  2 

/   (0.42+0.5  cos  -x  +  0.08  cos  — x)   dx 

-0  5                 a              a 
ENBW  =  -Olf^ 2  - 

[   /   (0.42  +0.5  cos-x  +  0.08  cos  — x)  dx] 
-0.5  a  a 

0.2646    .. 
"  0.1764  "  1,b 

0.5  ? 

COH.GAIN  =    /   (0.42  +  0.5  cos  -s  +  0.08  cos  — x)  dx  =  0.46 
-0.5  a  a 


The  window  and  its  transform  are  shown  in  Fig.  3.7(a)  and 

(b)  respectively.   The  main  lobe  width  is  3.0,  three  times  that 

of  the  rectangular  window,  and  its  first  side  lobe  level  is 

41  dB  down  from  the  main  lobe.   The  side  lobe  fall  off  rate  is 

18  dB/octave. 
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Figure  3 . 7 
(a)  Blackman  Window,  (b)  Spectral  Window 
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G.   KAISER  WINDOW  [7] 

The  Kaiser  window  is  also  called  a  Kaiser-Bessel  window 
and  its  mathematical  expression  is  given  as: 


I  CS/I.O-C-)2] 


w(k)  =  -2 — j  [g]  ° —  lxl^         (3.22) 


where  I  (x)  =  zeroth  order  modified  Bessel-function,  and  is 
o 

defined  as : 

,x.n   2 

00  (  ) 

I  (x)  =  1  +    EC  -V-  ]  (3.23) 

O  n  1 

n=l    n 


The  Fourier  transform  of  this  window  is  given  as 


9           sinhES  v4-(^-)2] 
W(f)    =    T    fflS      ra f<fa  (3. 24) 

1  K&  )  j Z ~ 


where  fa  =  $  /tt  . 

From  (3.24)  if  f=fa,  the  term  inside  the  square  root 

2 
becomes  zero  and  W(fa)  =  =? — r-o-\-      If  f>fa,  the  term  becomes 

V3) 


i  /<£>2-i, 


/(■jp— )    -I?    and   using   the    identity    sinh    ix    =    i    sin   x,    (3.24) 
becomes : 


9        sinCS  Z(^-)2-l] 
W(f)    =    t-Pq-y  ra   f >fa  •  (3.25) 

s/7X7-i 
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The  usual  range  of  3  is  between  4  and  9 ,  corresponding  to 
a  range  of  relative  peak  heights  from  3.1%  down  to  0.047%. 
As  3  increases,  the  side  lobes  decrease  in  amplitude  while 
the  main  lobe  widths  increase.   By  varying  3  the  proportion- 
ality between  the  width  of  the  main  lobe  and  the  height  of 
the  side  lobe  can  be  adjusted  as  desired. 

When  it  is  given  3=6  and  a=0.5  then  the  parameters  ENBW 
and  COH.GAIN  are  given  as: 


0.5   I  [6  A-C^)  ]  9 

/    (  °    j  u°rs     )2  dx 

ENBW  =  =°i* -^— =  °'1^7A    =  1.5 


K~~^^ 


0.25028 


0.5  IA[6A-(?p-)'"]      0 
v        c         °       0-5      jt2 

c    '   — rnn dx] 

-0.5       o 


0.5  ^L^1-^^ 

COH.GAIN  =     /    dx  =  0.5 

-0.5        I  [6] 
o 


Fig.  3.8(a)  and  (b)  shows  the  window  and  its  transform 
respectively.   The  main  lobe  width  is  2.0,  twice  that  of 
the  rectangular  window,  and  the  first  side  lobe  level  is 
43  dB  down  from  the  main  lobe.   The  side  lobe  fall  off  rate 
is  about  6  dB/octave. 

For  3=6  and  3=9  the  window  resembles  to  Hamming  and 
Blackman  window  respectively. 
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Figure  3 . 8 
(a)  Kaiser  Window,  (b)  Spectral  Window 
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TABLE  II I. 1 
ONE-DIMENSIONAL  WINDOW  AND  FIGURES  OF  MERIT 


WINDOW  ENBW     COH.GAIN     a±        a2     b 


dB  dB     Hz    dB/octave 

Rectangular  1.0  1.0  -13  -35  0.81  -6 

Barlett  1.33  0.5  -26  -58  1.63  -12 

Cosine  1.23  0.64  -23  -65  1.35  -12 

Hann  1.5  0.5  -32-88  1.87  -18 

Hamming  1.36  0.54  -43  -53  1.91  -6 

Blackman  1.57  0.46  -58  -94  2.82  -18 

Kaiser  1.5  0.49  -43  -66  2.2  -6 
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IV.   TRANSFORMATION  FROM  ONE-DIMENSION 
TO  TWO-DIMENSION  [4,  8,  9] 


The  two-dimensional  Fourier  transform,  two-dimensional 
window  functions,  and  window  figures  of  merit  are  discussed 
in  this  chapter. 

A.   TWO-DIMENSIONAL  FOURIER  TRANSFORM  [8] 

The  Fourier  transform  of  f(x,y),  a  function  of  two  inde- 
pendent variables  x  and  y,  is  .represented  by  F  {f(x,y)}  or 
F(u,v),  as  is  defined  as: 

oo    00 

F{f(x,y)}  =  F(u,v)  =   /   /   f(x,y)e"j(ux+vy)dxdy     (4.1) 


—  00   —  oo 


The  transform  is  itself  a  function  of  two  independent 
variables  u  and  v,  which  are  generally  referred  to  as  fre- 
quencies . 

Similarly  the  inverse  Fourier  transform  of  a  function 
F(u,v)  is  represented  by  F   (F(u,v)}  or  f(x,y)  and  is 
defined  as: 


00       00 


F_1{F(u,v)}  =  f(x,y)  =  /    /   F(u,v)ej(ux+vy)dudv    (4.2) 


_oo    — oo 


As  mathematical  operations,  the  transform  and  its  inverse  are 
very  similar  differing  only  in  the  sign  of  the  exponent 
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appearing  in  the  integral. 

In  order  that  the  integral  exist,  the  constraints  on 
the  function  f(x,y)  are  given  as: 

1.  f(x,y)  must  be  absolutely  integrable  over  the  infinite 
x,y  plane. 

2.  f(x,y)  must  have  only  a  finite  numbers  of  discon- 
tinuities and  a  finite  number  of  maxima  and  minima 
in  any  finite  extent. 

3.  f(x,y)  must  have  no  infinite  discontinuites  . 

The  basic  definition  of  the  Fourier  transform  pair  leads 
to  a  rich  mathematical  structure  associated  with  the  trans- 
form operation.   Considered  below,  a  few  of  the  basic  mathe- 
matical properties  that  will  find  use  in  later  discussion. 
These  properties  are  presented  as  mathematical  theorems, 
followed  by  a  brief  statement  of  their  physical  significance. 

1.   Linearity  theorem: 


F  (af+Sg>  =  aF  {f}  +  3  F  (g)  (4.3) 


The  transform  of  a  sum  of  two  functions  is  simply  the  sum 
of  their  individual  transforms. 
2.   Sealing  theorem: 


F  {g(ax,by)}=  -r^U  G  (£,  g)  (4.4) 


A  stretching  of  the  coordinates  in  the  space  domain  (x,y) 
results  in  a  construction  of  coordinates  in  the  frequency 
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domain  (u,v),  plus  a  change  in  the  overall  amplitude  of  the 
spectrum. 

3.  Shifting  theorem: 

t?  tef  ^w         Tf         \    -i(ua+vb)  ,.  _  N 

F  {f(x-a,  y-b)}  =  F(u,v)e  J  (4.5) 

A  translation  of  a  function  in  the  space  domain  intro- 
duces a  linear  phase-shift  in  the  frequency  domain. 

4.  Parseval's  theorem: 


00  00  00  oo 

/        /       |f(x,y)         dxdy    =      f        f       |F(u,v)|       dudv  (4.6) 

—  00  —00  —00  —00 


This  theorem  is  interpretable  as  a  statement  of  conserva- 
tion of  energy. 

5.  Convolution  theorem: 

00       00 

F  {   /    /  f (z,n)g(x-z,y-n)dzdn>  =  F(u,v)  G(u,v)    (4.7) 

—  00     —00 

The  convolution  of  two  functions  in  the  space  domain  is 
entirely  equivalent  to  the  simpler  operation  of  multiplying 
their  individual  transforms. 

6.  Autocorrelation  theorem: 


00       00  . 

F  {   /    /   f(z,n)f"(x-z,y-n)dzdn}    |F(u,v)|     < 4 • 8 > 

—  00     —00 
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similarly: 


00  CO  . 

F{|f(x,y)|     }    -      f        f      F(z,n)F"(z+u,n+v)dudv  (4.9) 


—  00  —  oo 


These  theorems  are  special  cases  of  the  convolution 
theorem. 

7.   Fourier  identity: 

At  each  point  of  continuity  of  f(x,y): 


F."1  xe  (f(x,y)}  =  F  F"1{f(x,y)}  =  f(x,y)        (4.10) 


The  successfive  transformation  and  inverse  transforma- 
tion of  a  function  yields  that  function  again,  except  at 
point  of  discontinuity. 

1.   Separable  and  Circularly  Symmetric  Functions  [9] 
A  function  of  two  independent  variable  is  called 
separable  with  respect  to  a  specific  coordinate  system  if  it 
can  be  written  as  a  product  of  two  functions. 

A  function  f(x,y)  is  separable  in  the  rectangular 
coordinate  x,y  if: 


f(x,y)  =  f  (x)  f2(y)  (4.11) 


A  function  f(r,90  is  separable  in  the  polar  coordinate 
r,9  if: 


f(r,&)  =  f  (r)  f2(90  (4.12) 
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Separable  functions  are  more  convenient  to  deal  with 
than  more  general  functions .   Separability  allows  complicated 
two  dimensional  manipulations  to  be  reduced  to  simpler  one 
dimensional  manipulations. 

A  function  separable  in  rectangular  coordinates  has 
the  simple  property  that  its  Fourier  transform  is  obtained  as 
a  product  of  two  one  dimensional  transforms 


OO       00 


F  (f(x,y)}  =   /    /   f(x,y)e~j(xu+vy)dxdy 


—  OO       00 


/   f1(x)e":)UXdx   /  f2(y)e"]vydy 

—  00  —CO 


F  {f  (x)}  F  (f2(y)}  (4.13) 


Thus  the  transform  of  f(x,y)  is  itself  separable  into 
a  product  of  two  factors.   The  process  simplifies  to  a  familiar 
one  dimension  manipulation  with  a  function  of  u,  and  a  function 
of  v. 

Functions  separable  in  polar  coordinates  are  not  so 
easily  handled  as  those  separable  in  rectangular  coordinates. 
But  it  is  generally  possible  to  demonstrate  that  two  dimen- 
sional manipulations  can  be  performed  by  means  of  a  series  of 
one  dimensional  manipulations.   The  simplest  class  of  functions 
separable  in  polar  coordinates  are  those  functions  which  ex- 
hibit circular  symetry. 
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A  function  is  said  to  be  a  circularly  symmetric  if  it 
can  be  written  as  a  function  of  radius  r,  independent  of  the 
variable  9 . 

f(r,9)  =  f(r)  (4.14) 

Such  functions  are  very  important,  since  a  window  function  can 
be  extended  into  a  two  dimensions  if  it  is  given  a  circularly 
symmetric  form.   Accordingly,  special  attention  should  be 
paid  to  these  functions. 

The  Fourier  transform  of  a  function  in  rectangular 
coordinates  is  given  by: 


OO       00 


F(u,v)  =   /    /   f(x,y)e~j(ux+vy)dxdy  (4.15) 


—  00     —00 


If  both  x,y  and  u,v  coordinates  are  transformed  into 
a  polar  coordinate  system  then: 


/  2   2 

r  =  /  x  +y  x=rcos| 


9  =    tan_1(-)  y=rsin9 

y 

and       

/  2^2  u  =  pcos4> 

o  =  /  u  +v 


$  =  tan   (— )  v  =  psin<J> 


Substituting  the  polar  coordinates  into  (4.15)  and 
using  the  transform  identity: 
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dxdy    =    rd9dr 

ux+vy   =    pcos<J)r   cos<f>  +  psin<{)r    sin9    =      r   cos(9-4>) 
The    Fourier    transform   in    (4.15)    becomes: 


F(p,<|>)    =    /         /    r  f(r,9)e~]pr   cos(9-4"    de    dr  (4.16) 

0  0 


If  the  Bessel  function  identity 


J  CprO  =  i-  T   e^pr  c°s(9-0)  dQ 
O        2tt   q 


is  used  to  simplify  (4.16),  then  the  transform  becomes: 


F(p)  =  2tt   /  r  f(r)  J  (pr)  dr  (4.18) 

Q 


Thus  the  Fourier  transform  of  a  circularly  symmetric 
function  is  itself  circularly  symmetric,  and  can  be  found  by 
performing  a  one  dimensional  manipulation  as  per  (4.18).   The 
expression  in  (4.18)  is  known  as  the  Fourier-Bessel  transform 
or  as  the  Hankel  transform. 

By  identical  arguments ,  the  inverse  Fourier  transform 
of  a  circularly  symmetric  function  F(p)  can  be  expressed  as: 


f(r)  =  2tt   /   pF(p)  J  (rp)  dp  (4.19) 

0         ° 
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There  is  no  difference  between  the  transform  and  its  inverse 
operation  in  the  case  of  Hankel  transform.   Using  the  notation 
H  {  }  to  represent  the  Hankel  transform  operation: 


H  H"1  (f(r)}  =  H  H  (f(r)}  =  f(r)  (4.20) 


From  the  Hankel  transform  pair  (4.18)  and  (4.19),  two 
dimensional  Fourier  transform  may  be  found  from  one  dimen- 
sional mainpulations . 

Huang  [4]  has  shown  that  a  good  two  dimensional  window 
function  can  be  derived  from  a  good  one  dimensional  function, 
as  long  as  it  is  circularly  symmetric.   This  is  discussed  in 
the  next  section. 

The  Hankel  transform  in  (4.18)  will  be  used  to  cal- 
culate the  Fourier  transform  of  a  two  dimensional  window. 

B.   TWO  DIMENSIONAL  WINDOW  FUNCTION  [4] 

If  W?(u,v)  and  w?(x,y)  is  a  two  dimensional  Fourier 
transform  pair,  then  its  inverse  transform  relation  is  given 
as  : 


w  (x,y)  =   /    /   W2(u,v)e:(uX+vy)dvdu  (4.21) 

-co     —00 


For  y=0,  (4.21)  becomes: 


00       CO 


w2(x,o)  =   /  {  /  W2(u,v)  dv}e]UX  du  (4.22) 


-co    —co 
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If  w?(x,y)  is  given  as  a  circularly  symmetric  function, 

/  2   2 
then  w?(x,y)  =  w(  /  x  +y  ).   For  circularly  symmetric  function 


w2(x,o)  =  w(x)  (4.23) 


And  its  inverse  Fourier  transform  is  given  as: 

w(x)  =   /  W(u)e:ux  du  (4.24) 


—  CO 


From  (4.22)  and  (4.23)  it  is  evident  that 


W(u)  =   /  W  (u,v)dv  (4.25) 

—  00 


The  expression  in  (4.25)  is  called  a  profile  function, 
and  similarly  the  profile  function  in  v  is  expressed  as: 


W(v)  =   /  W2(u,v)du  (4.26) 


Now  consider  the  following  two  equations: 


1     for   u  _>  o 
H(u)  =   {  (4.27) 

0         u  <  o 


and 


1   for  u  _>  o  for  all  v 
H  (u,v)  =  {  (4.28) 

0        u  <  o  for  all  v 
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The  expressions  in  ( M- .  2 7 )  and  (4.28)  are  one-dimensional  and 
two  dimensional  step  function  respectively. 

If  W?(u,v)  is  a  two  dimensional  window's  Fourier  trans- 
form, then  its  convolution  with  the  step  function  is  given 
as: 


00      oo 

W2(u,v)8  H  (u,v)  =   /    /   W2(t,s)H(u-t,  v-s)dtds 

OO       00 


U        °° 

/   {   /  W2(t,s)ds}dt  (4.29) 

o    -°° 


where  t,s  are  dummy  variables. 

Using  the  profile-function  identity  given  in  (4.25)  and 
(4.26),  it  can  be  shown  that: 


W(t)  =   /  W  (t,s)ds  (4.30) 

—  oo 


Equation  (4.29)  becomes 


u 
W2(u,v)  0  H2(u,v)  =   /  W(t)dt  (4.31) 

o 


Expanding  the  right  hand  side  of  (4.31)  into  a  convolution  form 


u  °° 

/  W(t)dt  =   /  W(t)  H(u-t)dt 


=  W(u)  0  H(u)  (4.32) 
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Then  substituting  (4.32)  to  (4.31): 


W  (u,v)  0  H2(u,v)  =  W(u)  0  H(u)  (4.33) 


Huang  states  that  if  a  two  dimensional  function  is  given 
in  circularly  symmetric  form,  then  its  Fourier  transform 
satisfies  the  relation  given  in  (4.33). 

To  show  the  benefit  of  the  results  in  (4.33)  consider  an 
example  of  designing  a  low  pass  filter. 
EXAMPLE : 

The  ideal  Fourier  transform  of  one  dimensional  Low  Pass 
Filter  is  given  as: 


1    for  |u|  _<  A 
F(u)  =  {  (4.34) 

0        lul  >   A 


The  impulse  response  f(x)  of  the  filter  is  equal  to  its 
inverse  Fourier  transform,  which  has  an  infinite  duration. 
In  reality,  the  impulse  response  has  to  be  finite.   By  using 
a  window  function,  the  infinite  duration  of  the  impulse  res- 
ponse is  truncated  to  become  a  finite  impulse  response  as 
given  as : 

g(x)  =  f(x)  •  w(x)  (4.35) 


where 


w(x)  =  0       for  |x|>'a  (4.36) 
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The  windowed  frequency  response  of  the  FIR  low  pass  filter 
is  now  equal  to  the  convolution  F(u)  and  W(u) : 

G(u)  =  F(u)  0  W(u)  (4.37) 

The  ideal  Fourier  transform  of  a  two  dimensional  LPF  is 
given  as  : 


1       |u2+v2|_<A2 
F.(u,v)  =  {         o   o  '   o  (4.38) 


0       |u2+v2|>A2 


The  infinite  impulse  response  f_(x,y)  of  this  filter  has 
to  be  truncated  using  a  window  function  w  (x,y).   The  finite 


impulse  response  filter  is  given  as: 


g  (x,y)  =  f  (x,y)  •  w0(x,y)  (4.39) 


The  windowed  frequency  response  of  this  filter  is  given  as 


G2(u,v)  =  F2(u,v)  8  W2(u,v)  (4.40) 


A  good  window  should  have  the  width  of  its  frequency  response 
much  smaller  than  the  width  of  the  filter's  frequency  response 


72 2 
W(u)  and  W?(u,v)  should  be  zero  for  values  of  /u  +  v  much 

smaller  than  A. 
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Near  a  discontinuity  in  the  filter's  frequency  response 
F(u),  i.e.,  u  =  +A,  G(u)  is  essentially  a  convolution  of 
W(u)  and  a  step  function: 

G(u)  =  H(u)  8  W(u)  (4.41) 

Similarly  near  a  discontinuity  in  the  filter's  frequency 

*  2   2     2 

response  F  (u,v),  i.e.,  u  +v   =  A  ,  G2(u,v)  is  essentially  a 

convolution  of  W~(u,v)  and  a  step  function: 


G2(u,v)  =  H  (u,v)  0  W2(u,v)  (4.42) 


From  Huang's  theorm  given  in  (4.33),  G_(u,v)=G(u)  if  only 

2"   2 
w  (x,y)  =  w(  /  x  +y  ) . 

This  means  that  if  w(x)  is  a  good  window,  then,  w  (x,y) 
is  also  a  good  two  dimensional  window. 

If  other  words,  a  good  two  dimensional  window  can  be 
derived  from  a  good  one  dimensional  window  by  extending 
it  into  a  circularly  symmetric  function. 

C.   TWO  DIMENSIONAL  WINDOW  FIGURE'S  OF  MERIT 

The  window's  figures  of  merit  found  in  the  one  dimen- 
sional case  are  now  applied  to  the  two  dimensional  case.   The 
same  figures  of  merit  that  may  be  used  for  comparison  among 
windows,  will  be  re-examined  to  determine  if  there  should 
be  any  changes  for  the  two  dimensional  case. 
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1.   Equivalent  Noise  Band  Width  (ENBW) 

The  ENBW  in  one  dimensional  case  is  defined  as  the 
width  of  a  rectangular  filter  with  the  same  peak  power  gain 
that  passes  the  same  noise  power. 

A  rectangular  filter  in  two-dimension  represents  a 
cross-section  of  a  cylindrical  filter.   The  noise  power 
passed  by  the  area  under  the  rectangular  filter,  becomes 
noise  power  passed  by  the  volume  under  the  cylindrical  filter 

If  the  peak  power  gain  is  the  same,  then  the  volume 
divided  by  the  height,  results  in  the  area  covered  by  the 
base  of  the  cylinder,  which  is  a  circle. 

Two-dimensional  ENBW  is  defined  as  the  circumference 
of  the  horizontal  cross-section  of  a  cylindrical  filter, 
that  passes  the ■ same  noise  power. 

The  vertical  cross-section  of  a  two  dimensional 
window  is  equal  to  a  one  dimensional  window.   Consequently 
the  width  of  the  equivalent  rectangular  window  is  equal  to 
the  diameter  of  the  equivalent  cylindrical  window.   Then 
the  two  dimensional  ENBW  is  given  as: 

(ENBW)2D  =  tt(ENBW)1d  (4.42) 

For  the  two  dimensional  case,  the  definition  of  ENBW 
should  be  changed  from  the  width  of  a  rectangular  to  the 
circumference  of  a  cylinder,  and  the  value  can  be  calculated 
directly  from  the  one  dimensional  case. 
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2 .   Coherent  Gain 

Two  dimensional  coherent  gain  is  best  defined  using 
two  dimensional  sinusoidal  signal.   The  input  signal  is 
defined  as: 

f(x,y)  =  Av   ej(ux+vy)  +  n(x,y)  (4.43) 

x ,  y 

where  n(x,y)  is  two  dimensional  white  noise. 

The  output  signal  or  the  windowed  signal  is  given  as 

f  (x,y)  =  Aw(x,y)ej(ux+vy)+w(x,y)n(x,y)        (4.44) 
w 

The  two  dimensional  frequency  response  of  the  win- 
dowed signal  is  given  as: 

v    t         \  r         r      a  r         \    j(ux+vy)  -j(ux+vy),  , 

F  (x,y)  =  f         J      Aw(x,y)eJ      J    e  J      J    dxdy 
w 

-a  -a 


a    a 
=  A  /    /  w(x,y)dxdy.  (4.45) 

-a   -a 


The  amplitude  of  the  windowed  signal  frequency  response 
is  proportional  to  the  original  amplitude  A.   The  propor- 
tionality factor  is  the  double  integral  of  the  window  function. 

If  the  window  function  is  given  in  circularly  symme- 
tric form,  then  the  proportionality  factor  becomes: 

a    a  a   2tt 

/    /  w(x,y)dxdy  =   /    /   w(r)rd9dr 
-a   -a  o    o 
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a  2tt 

/  r  w(r)dr   /  dl 
o  o 


a 
=  2tt   /  r  w(r)dr  (4. 46) 

o 


The  proportionality  factor  for  a  rectangular  window 

2 
is  equal  to  Tra  .   The  proportionality  factor  for  the  other 

window  is  reduced  due  to  the  window  smoothly  going  to  zero 

near  the  boundaries . 

If  the  proportionality  factor  in  (4.46)  is  normalized 

2 
by  factor  iTa   (that  is,  divided  by  the  area  covered  by  the 

base  of  the  window) ,  then  the  proportionality  factor  becomes 

coherent  gain: 


2   a 
COH.GAIN  =  -4-  /  r  w(r)dr  (4.47) 


The  coherent  gain  given  in  (4.47)  above  is  used  as 
window's  figure  of  merit. 

3 .   Trade-off  Parameters 

Two  dimensional  trade-off  parameters  are  similar  in 
definition  with  the  ones  used  in  the  one  dimensional  case. 
The  trade-off  parameter  a, ,  a~  ,  b  and  d  are  derived  from  the 
frequency  response  of  a  window,  which  is  found  using  the 
Hankel  transform. 

The  Hankel  transform  allows  a  two  dimensional  trans- 
form to  be  found  from  a  one  dimensional  manipulation.   Con- 
sequently the  definitions  of  the  trade-off  parameters  doesn't 
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change.   The  graphical  interpretation  of  the  parameters  are 
shown-  in  Fig.  4.1, 
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Plot  Of  Kaiser  Window,  Used  To  Illustrate  The 
Definition  Of  Trade-Of f  Parameters  a, ,  a« ,  b  and  d 
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V.   TWO  DIMENSIONAL  WINDOW  FUNCTIONS 

The  two  dimensional  extension  of  the  previous  window 
functions  will  be  given  in  this  chapter.   The  window's  fre- 
quency response  is  found  using  Hankel  transform.   Unfor- 
tunately it  is  hard  to  present  the  Hankel  transform  of  a 
window  in  a  closed  form.   Only  two  transforms  are  presented 
in  closed  form,  they  are  Rectangular  and  Barlett  windows, 
the  rest  are  presented  in  their  original  integral  form. 

Without  loss  of  generality,  the  window  function  is  given 
as  w(r) ,  a  two  dimensional  circularly  symmetric  function, 
restricted  to  the  interval  |r|  _<  a ,  while  the  window  is  set 
to  unity  at  the  origin. 


00 

w(o)  =   /  f  W(f)  df  =  1  (5.1) 

—  00 


w(r)  =  o  for   |r|  >  a  (5.2) 

The  frequency  response  of  this  window  is  given  as: 


W(f)  =   /  r  w(r)  J  (fr)  dr  (5.3) 

o  ° 


The  frequency  response  consists  of  a  main  lobe  and  side 
lobes  which  posses  a  large  amplitude  ratio  in  a  linear  scale 
The  frequency  response  plot  is  given  in  log-frequency  scale 
versus  log-normalized  amplitude  (W(f)/W(0)). 
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A.   RECTANGULAR  WINDOW  [2] 

The  two  dimensional  version  of  the  circularly  symmetric 
rectangular  window  is  given  as  follows: 


1      for  |  r  |  _<  a 
w(r)  =  {  (5.4) 

0     for  Irl  >   a 


This  window  is  cylindrical,  with  the  height  equal  to  one  and 
the  diameter  of  the  base  equal  to  2a.   The  vertical  cross- 
section  of  this  window  is  equal  to  a  rectangular  window. 

The  frequency  response  of  this  window  is  given  by  the 
Hankel  transform  as : 


a 
W(f)  =   /  r  J  (fr)  dr  (5.5a) 

o 


The  recurrence  formulae  is  used  to  solve  the  integral 
in  (5.5a),  where:   [10] 


-r—    {x  J  ( x )  }   =xJ   ,  ( x ) 
dx     n  n-1 


—-  {x  J. (x)}   =  x  J  (x) 

dx     1  o 


Then  the  integral  form  in  (5.5a)  reduces  to: 

1   a 
W(f)  =  -±-     I    (fr)J  (fr)  d(fr) 

f2   o      ° 

=  -%  /  d(fr)J  (fr)  =  {-^-   J  (fr)  | 
f   o  f .         o 
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2  J,(af) 

=  a   — —5 .  (5  .  5b) 

J  (af) 
The  window  frequency  response  { — -, — ^r-}  has  a  form  similar  to 

s  inx 
( )  form.   The  value  of  this  function  at  the  origin, 

f=0,  is  one.   The  envelope  of  the  window  frequency  response 

crosses  the  frequency  axis,  whenever  J,(af)=0.   The  first 

zero-crossing  occurs  at  the  value  of  af=3.8.   If  the  value 

of  a  is  increased  (meaning  the  window  impulse  response  is 

lengthened) ,  the  zero-crossing  occurs  at  smaller  value  of 

f  (meaning  the  window  frequency  response  is  narrower) .   If 

the  value  of  a  is  decreased  (meaning  the  window  impulse 

response  is  shortened) ,  the  zero-crossing  occurs  at  larger 

value  of  f  (meaning  the  window  frequency  response  is  wider) . 

At  the  value  of  a=0.5,  the  first  zero-crossing  occurs  at 
the  value  of  f=7.6. 

The  ENBW  of  this  window  with  the  value  of  a  is  given  as 
a=0.5: 


ENBW  =  tt(1)  =  3.14 


The  COH.GAIN  is  given  as: 

COH.GAIN  =  — — rr      f      r  dr  =   ^rl   =1.0 

7T(  0  .  5  )      O 


The  window  and  its  frequency  response  are  given  in 
Figure  5.1(a)  and  (b)  respectively.   The  main  lobe  width 
(the  width  between  the  origin  and  the  first  zero-crossing)  is 
7.6.   The  first  side  lobe  level  in  18  dB  down  from  the  main 
lobe  peak  and  the  side  lobe  fall-off  rate  is  9  dB/octabe. 
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l  Figure  5 , 1 

(a)  2-D  Rectangular  Window,  (b)  2-D  Spectral  Window 
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B.   BARTLETT  WINDOW 

The  two  dimensional  circularly  zymetric  Bartlett  window 
is  given  as : 

1    1  g  1  |  r  |  _< 

w(r)  =  {      a  (5.7) 

0  elsewhere 

The  shape  of  this  window  function  is  similar  to  a  cone, 
with  the  height  equal  to  one,  and  the  base  in  a  circle  with 
a  radius  r=a.  The  vertical  cross-section  of  this  window  is 
equal  to  a  one-dimensional  Bartlett  window. 

The  frequency  response  of  this  window  is  given  by  the 
Hankel  transform  as  follows: 

a  I     I 

W(f)    =      /      r(l    -   -LEI)    J    (fr)    dr 

a  o 

o 


a  a  2 

/      r   J    (fr)dr  -      /  —  J    (fr)    dr 

o  a       o 

o  o 


J    (af)         x      a        d    (fr)J    (fr) 

=    a   — a      f   r   ~2 

o  f 


J,(af)              Jn(af)         ,       a    (fr)J    (fr) 
a   ^—  -    a   -V-  +    |      /   f2 dr 


n       a    (fr)2J, (fr) 

|     /   nji d(fr) 

o  f   r 


,       a   d(fr)2J9(fr) 


a      J  PU 

o  f  r 
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i  J0(af)    -  i     /    (fr)2   J9(fr)    d    C-jL) 
rr      2  a  2  ^4 

o  f  r 


J9(af)  a 

+   -At     /   J«(fr)    d(fr) 


.2  .3      J    "2 

af        o 


J9(af)  ,  a 

-^-75 —  +   -^y    {      /   J    (fr)    d(fr)    -    2J,(af)} 
r  af  o 


1  1      a 

{af   Jn(af)    -    2    J,  (af)}    +    -Ar  /   J    (fr)    d(fr) 


3     ^    -2—         -    TB1,J  .3    J    uo 

af  af      o 


l  2  a 

{-af   J    (af)}    +       a     ,      /   J    (fr)    d(fr) 


,-3  O  /         r^3  o 

af  (af)         o 


W(f)    =    a2{ — i-,y  J    (af)    +    — i-_.      /   J    (fr)    d(fr)}         (5.8) 
(af)  (af)         o 


If  p=af  is  substituted  into  (5.8)  then  the  end  result  is: 


W(p  )  =  a2{p"3  /  J  (t)dt  -  p'2J  (p)}  (5.9) 

a  o  o 

o 


where  t  is  a  dummy  variable  of  integration. 

The  window  and  its  frequency  response  is  given  in  Figure 
5.2(a)  and  (b)  respectively,  where  a=0.5.   The  ENBW  and 
COH.GAIN  parameters  are  given  as: 

ENBW  =  tt(1.33)  =  4.18 

COH.GAIN    =    — — *•      /      r(lr    j^|    )    dr    =    ~||   =    °-34 

M(0.5)2      o  0>5  °'78 
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Figure  5 . 2 
(a)  2-D  Barlett  Window,  (b)  2-D  Spectral  Window 
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The  main  lobe  width  is  about  11.5,   The  first  side  lobe 
level  is  34  dB  down  from  the  main  lobe  level,  and  the  side 

lobe  fall  off  rate  is  12  dB/octave. 

P 
The  term   /  J  (t)dt  in  (5.9)  is  hardly  a  closed  form, 

o   ° 
but  it  is  tabulated.   Then  the  rectangular  and  relatively 

the  Bartlett  windows  are  two-dimensional  windows  that  have 

a  closed-form  frequency  response. 

C.   COSINE  WINDOW 

The  third  two-dimensional  window  to  be  discussed  is  a 
cosine  window.   The  two  dimensional  circularly  symmetric 
window  is  given  as : 

T  c  I   I 

cos  ~—  r  for   |r|  <  a 

w(r)  =  {      Za  "  (5.10) 

o  elsewhere 

The  shape  of  this  window  is  similar  to  a  half  sphere  and 
the  height  at  r=o  is  one.   The  base  of  the  window  is  a  circle 
with  radius  r=a. 

The  frequency  response  of  the  window  is  given  as  a 
Hankel  transform  with  no  closed  form. 


a 

w(f)  =   /  r  cos  #-  r  J  (fr)  dr  (5.11) 

2a    o 
o 


The  effort  to  derive  a  closed  form  is  given  in  the  Appendix 
If  a=0.5,  then  the  window  and  its  frequency  response  are 
given  in  Figure  5.3(a)  and  (b)  respectively.   The  parameters 
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Figure  5 . 3 

(a)  2-D  Version  Of  Cosine  Window  (cross-section), 
(b)  2-D  Spectral  Window 
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ENBW  and  COH.GAIN  are  given  as: 
ENBW  =  tt(1.33)  =  3.86 


COH.GAIN  =  — — o   /   r  cos  Trr  dr  =  ^44  =  0.46 

tt(  0  .  5  )    o 


The  main  lobe  width  is  10.5  while  the  first  side  lobe 
level  is  2  6  dB  down  from  the  main  lobe  and  the  side  lobe 
fall  off  rate  is  12  dB/octave. 

D.   HANN  WINDOW 

The  two  dimensional  circularly  symmetric  Hann  window 
is  given  as: 

cos  yr- r  =  0.5  +  0.5  cos  —  r      |r|  <a 
w(r)  =  {      Za  a  (5.12) 

o  elsewhere 

The  vertical  cross-section  of  the  window  is  equal  to  a 
one  dimensional  Hann  window. 

The  frequency  response  is  given  by  the  Hankel  transform 
as  : 

o 

w(f)  =   /  r(0.5+0.5  cos  -  r)  J  (fr)  dr 

a      o 
a 

J.(af)        a 

=  0.5  a— ±-s +0.5   /   cos  -  r  J  (fr)  dr     (5.13) 

f  a    o 

o 
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The  first  term  in  (5.13)  has  a  closed  form  as  in  the  case 
of  a  Rectangular  window.   The  second  term  has  no  closed  form 
as  in  the  case  of  the  cosine  window. 

For  a  window  length  of  a=0.5,  the  window  and  its  frequency 
response  are  given  in  Figure  5.4(a)  and  (b)  respectively. 
The  ENBW  and  COH.GAIN  parameters  are  given  as: 

ENBW  =  it  (1.5)  =  4.7 


0  .  5 
COH.GAIN    =    — *■      /   r   cos2    Trr   dr    =    ~|i   =    0.2  9 

—  /  n       c  \  ^  ~  U./O 

tt(  0  .  5)    o 


The  main  lobe  width  is  13.7  and  the  first  side  lobe  level 
is  34  dB  down  from  the  main  lobe  level.   The  fall-off  rate 
of  the  side  lobes  is  18  dB/octave. 

E.   HAMMING  WINDOW 

The  two  dimensional  circularly  symmetric  Hamming  window 
is  given  as : 

0.54  +  0.4  6  cos  -r  |r|  <  a 

w(r)  =  {  a  (5.14) 

o  elsewhere 

The  shape  of  this  window  is  a  smooth  cone  sitting  on  a 
cylindrical  base.   The  vertical  cross-section  of  the  window 
is  equal  to  the  one  dimensional  case. 

The  frequency  spectrum  of  the  window  is  given  by  the 
Hankel  transform  as: 
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Figure  5  .  M- 
(a)  Von  Hann  Window  in  2-D,  (b)  Spectral  Window 
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a 
w(f)  =   /  r(0.54  +  O.146  cos  -r)  J  (fr)  dr 
o 


J  (af)         a 

=  0.54  a  -^-= +  0 .  M-6  /  r  cos  -r  J  (fr)  dr     (5.15) 

f  a   o 

o 


The  first  term  in  (5.15)  is  the  familiar  frequency  re- 
sponse of  the  rectangular  window  and  the  second  term  is  the 
response  of  a  cosine  window. 

If  the  length  of  the  window  is  a=0.5,  then  the  window 
and  its  frequency  response  are  given  in  Figure  5.5(a)  and  (b) 
respectively.   The  ENBW  and  COH.GAIN  parameters  are  given  as: 

ENBW  =  tt(1.36)  =  4.3 


COH.GAIN  =  — — *■   /   r(0.54  +  0.4  6  cos  7Sr)dr=^4l=0  •  36 

*  n      c  \  £.  U.D        U./O 

tt(  0  .  5  )    o 


The  main  lobe  width  is  equal  to  13.5  and  the  first  side 
lobe  level  is  4  7  dB  down  from  the  main  lobe  level.   The  side 
lobe  fall  off  rate  is  6  dB/octave. 

F.   BLACKMAN  WINDOW 

The  two  dimensional  circularly  symmetric  Blackman  window 
is  given  as : 

0.42  +  0.5  cos—  +  0.08  cos—  r     |r|   :  a 
w(r)  =  {  a  a  (5.16) 

o  elsewhere 
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Figure  5 . 5 
(a)  Hamming  Window  in  2-D,  (b)  Spectral  Window 
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The  vertical  cross-section  of  this  window  is  equal  to 
the  one-dimensional  Blackman  window. 

The  window  frequency  response  is  given  by  the  Hankel 
transform: 


a 
W(f)  =  f   r{0.42  +  0.5  cos-  r  +  0.08  cos—  r}J  (fr)dr 
o  a  a     o 


J  (af)      a 

=  0.4  2  a-i-s +0.5  /  r  cos—  J  (fr)dr 

x  a   o 

o 


a        2tt 

+  0.08   /  r  cos  — r  J  (fr)dr      (5.17) 

a    ° 


If  the  extent  of  the  window  given  as  a=0.5,  the  window 
and  its  frequency  response  are  given  in  Figure  5.6(a)  and  (b) 
respectively.   The  ENBW  and  COH.GAIN  parameters  are  given  as: 

ENBW    =    tt(1.57)    =    4.9 


9  0.5 

COH.GAIN    =    ^V-      /      r    (0.42  +  0.5    cost^+0.0  8    cos?f^rr)dr 

(0.5)    it       o 


17 

=    0.22 


0.78 


The  main  lobe  width  is  equal  18.5,  which  is  almost  triple 
of  that  of  the  rectangular  window.   The  first  side  lobe  level 
is  60  dB  down  from  the  main  lobe  and  the  side  lobe  fall  off 
rate  is  9  dB/octave. 
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Figure  5 . 6 
(a)  2-D  Blackman  Window,  (b)  Spectral  Window 
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G.   KAISER  WINDOW 

Two  dimensional  circularly  symmetric  Kaiser  window  is 
given  as : 


£- 


i  [3  /l-c-^r] 


o       a 

w(r)  =  {       o^"1  (5.18) 

o  elsewhere 

The  vertical  cross-section  of  the  window  is  equal  to  a 
one  dimensional  Kaiser  window. 

The  frequency  spectrum  of  the  window  is  given  by  the 
Hankel  transform: 

I  [S/i_(M)2 


jO 


W(f)  =   /   r  — 1  |-  j JQ(fr)dr  (5.19) 


which  cannot  be  reduced  to  a  closed  form. 

If  a=0.5  and  8=5  then  the  window  and  its  frequency  res- 
ponse are  given  in  Figure  5.7(a)  and  (b) .   The  ENBW  and 
COH.GAIN  parameter  are  given  as: 

ENBW  =  tt(1.5)  =  4-  .71 


COH.GAIN  =       2    /   r  T~m dr  =  ^-^  -  0.61 

tt(  0 .  5 )     o         o 


The  main  lobe  width  is  14.0  and  the  first  side  lobe  level 
is  4-4  dB  down  from  the  main  lobe  peak.   The  side  lobe  fall  off 
rate  is  6  dB/octave. 

In  this  case,  the  window's  characteristics  depend  strongly 
on  the  values  of  a  and  3.   With  3=6,  the  window  is  similar  to 
a  Hamming  window. 
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Figure  5 . 7 
(a)  2-D  Kaiser  Window,  (b)  Spectral  Window 
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TABLE  VI. 1 
Two-Dimensional  Window  Functions  and  Figures  of  Merit 


WINDOW      ENBW  COH.GAIN  a  a  b  d 

dB  dB  Hz  dB/octave 

1.  Retangular  3.14  1.0  -18  -38  6.6  -9 

2.  Bartlett    4.18  .34  -34  -58  11.0  -12 

3.  Cosine      3.86  .46  -26  -55  10.0  -12 

4.  Hann        4.71  .29  -34  -64  13.0  -18 

5.  Hamming     4.3  .36  -47  -51  13.0  -6 

6.  Blackman    4.9  .22  -60  -66  19.0  -9 

7.  Kaiser      4.7  .61  -44  -61  14.0  -6 
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VI.   CONCLUSION 

One  of  the  most  straightforward  applications  of  the  rec- 
tangular window  is  to  truncate  an  infinite  length  impulse 
response.   This  can  be  thought  as  multiplying  the  impulse 
response  by  a  window  function  that  is  unity  in  some  finite 
interval  between  ±a,  and  zero  elsewhere.   This  multiplica- 
tion is  equivalent  to  convolving  the  Fourier  transform  of 

s  inx 
the  impulse  response  with  a  function  of  the  form  .   This 

causes  ringing  in  the  original  frequency  response.   Increasing 

the  length  of  the  finite  impulse  response  (by  increasing  the 

length  of  the  rectangular  window)  decreases  the  duration  of 

the  ringing,  but  has  no  effect  on  the  peak  overshoot.   This 

effect  is  called  the  Gibbs  phenomenon. 

A  considerable  amount  of  effort  has  been  put  into  finding 

a  good  window  function,  i.e.,  window  functions  which  reduce 

the  peak  overshoot  of  the  Gibbs  oscillation.   It  is  generally 

agreed  that  the  frequency  response  of  a  good  window  should 

have  a  narrow  central  lobe  (with  a  large  amplitude)  and  side 

lobes  with  comparatively  small  amplitude.   There  are  many 

window  functions  readily  available  in  the  literature,  however 

most  window  functions  are  given  only  in  one  dimension.   The 

one  dimensional  window  functions  and  their  Fourier  transform 

are  given  in  detail  in  Chapter  III  of  this  thesis.   Figures 

of  merit  for  one  dimensional  window  are  powerful  measure  and 

are  used  to  compare  the  characteristic  of  the  various  windows. 
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It  has  been  the  objective  of  this  thesis  to  examine  two 
dimensional  window  functions.   The  two  dimensional  window 
functions  presented  here  have  been  a  direct  extension  of 
corresponding  one  dimensional  window  functions.   The  two 
dimensional  window  function  is  formed  from  a  one  dimensional 


2   2 
window  function,  using  the  relation  w„  (x,y)  =w( /x  +y  ).   This 

obviously  rotates  the  one  dimensional  function  about  the 

origin  in  the  x-y  plane.   However  it  does  not  simply  rotate 

the  frequency  response  of  the  one  dimensional  about  the  origin 

in  the  u,v  plane.   The  two  dimensional  frequency  response  of 

the  window  can  be  calculated  using  the  Hankel  transform.   The 

principal  advantage  of  this  approach  is  that  the  frequency 

transform  of  two  dimensional  window  can  be  found  using  a  one 

dimensional  manipulation. 

The  two  dimensional  window  figures  of  merit  are  a  direct 
extension  of  the  one  dimensional  case.   However  the  definition 
of  the  ENBW  is  changed  from  the  width  of  a  rectangular  filter 
in  the  one  dimensional  case  to  the  circumference  of  a  circular 
filter  in  the  two  dimensional  case.   Table  VI . 1  shows  the 
comparison  between  one  dimensional  and  two  dimensional  window 
figure  of  merit.   From  the  table  it  is  clear  that  the  two 
ENBW  parameters  are  related  by  the  equation  (M-.M-2). 

The  coherent  gain  decreases  slightly  except  for  the  Kaiser 
window  where  it  is  higher,  and  the  rectangular  window  where  it 
is  constant.   The  first  side  lobe  level,  which  is  shown  by  the 
parameter  a..  ,  is  lower  for  the  two  dimensional  windows,  however 
the  width  of  the  main  lobe  is  much  greater  in  each  case.   The 


85 


side  lobe  fall-off  rate  stays  the  same  except  for  the  Blackman 
window  where  it  is  slower  and  the  rectangular  window  where  it 
is  faster. 

In  conclusion,  while  one  dimensional  and  two  dimensional 
windows  can  be  simply  related  to  one  another  in  the  time 
domain,  no  direct,  general  parallels  exist  to  relate  their 
transforms  in  the  frequency  domain. 
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TABLE  VI. 1 

Comparison  Between  Two  Dimensional  And 
One  Dimensional  Window  Figures  of  Merit 


ENBW      COH.GAIN      a 


WINDOW 


2D     ID    2D    ID     2D    ID    2D     ID      2D    ID 


Rectangular  3. 14  1.0  1.0  1.0  -18  -13  6.6   0.81  -9  -6 

Bartlett  4.18  1.33  .34  .5  .34  -26  11.0  1.63  -12  -12 

Cosine  3.86  1.23  .46  0.64  -26  -23  10.0  1.35  -12  -12 

Hann  4.71   1.5  .29  0.5  -34  -32  13.0  1.87  -18  -18 

Hamming  4.3   1.36  .36  0.54  -47  -43  13.0  1.91  -6  -6 

Blackman  4.9   1.57  .22  .46  -60  -58  19.0  2.82  -9  -18 

Kaiser  4.7    1.5  .61  .49  -44  -43  14.0   2.2  -6  -6 
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APPENDIX 

The  effort  to  put  a  Hankel  transform  of  a  cosine  function 
in  a  closed  form  is  given  below.   The  cosine  function  arises 
from  Hankel  transform  of  a  cosine  window. 

The  Hankel  transform  of  a  cosine  function  is  given  below: 


w(f)  =  /  r  cos  ~—  r  J  (rf)  dr        |r|  <  a       (A.l) 

za    o  ii— 

o 


The  above  integral  form  has  no  known  closed  form,  as  proved 
by  the  following  discussion. 

Three  different  methods  are  used  to  bring  the  integral 
into  a  closed  form,  they  are:   straight  integration,  partial 
differential  and  table  of  integration. 

A.   STRAIGHT  INTEGRATION 

To  use  straight  integration  method  to  solve  the  integra- 
tion, the  variable  needs  to  be  changed  into  trigonometric 
variables . 

T      a    irr        2a 
Let :   9  =  o—   ,  r  =  —   6  . 
2a   '      it 

TT 

The  interval  of  the  integral  becomes,  9=o  and   =y  for  r=o 

and  r=a  respectively. 

Substitution  of  the  new  variables  into  Hankel  transform 

gives : 

X  2  */2 

w(f)  =  (^)    /   e  cos  e  J  (xe)  de  (A. 2) 

f  o 

o 
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where  A  =  2( — ).   This  new  integral  form  gives  more  pleasing 

appearance  than  the  previous  one. 

Substitute  the  Bessel  function  J  (A9)  inside  the  integral 

o  & 

with  Bessel  equation,  given  as  follows: 


J  "(x)  +  -  J  '  (x)  +  J  (x)  =  0  (A. 3) 

O        X   o  o 


Let  x  =  A9,  and  do  some  manipulation.   The  Bessel  equation 
becomes : 


J  (X9)  =  -  4y  {J  M(A9)  +  i  J  »(A9)}  (A. 4) 

o         ,  /    o         6   o 


The    integral    form   in   A. 2    becomes 


1     ff/2  1 

w(A)    =    ~      f    9    cos    9     {J    M(A9)    tij    '(A9)}df 

f2      o  °  9      ° 


1  tt/2  it/2 

~{9    cos  9    J    '(A9)  /    {(cos    9-9    sin    9)J    '(A9) 

cZ  °  «  O 

X  o  o 


-    cos    9    J    ? (A9)}d9} 
o 


,  tt/2 

•-4-    {      /      9    sin    9    J    '  (A9)    d9} 
f2 


,  tt/2         tt/2 

~  {9    sin    9    J    (A9)  -      /      (sin    9+9    cos    9)J    (A9)d9> 

r-Z  O  '  O 

f  o  o 


1  1     7r^2                                         1 

•^-{3-  J    (5-A)}  +    -4-      /  sin    9    J    (A9)    d9    +    4=r  w(A) 

ct-    I      o    I  ci                                o                           ,2 

f  f        o                                                  A 
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1    Tt/2 

(1    -  4r)wU)    =    ~£-  J    (Jx)    +   ^     /    sin    8    J    (A9)    df 

X2  2f2      °    2  f2      o 


The  end  result  is: 


,2         tt/2  ,,2       . 

w(A)    =    -V1^      /    sin    8    J    (A6)    d6    -    ^    :1;    J       (£x)       (A. 5) 
A2f2      o  °  2A2f2  °      2 


By  replacing  J  (A9)  inside  the  integral  as  before,  it  may 
be  continued  at  will,  but  the  term  with  the  integral  form 
cannot  be  eliminated.   The  conclusion  is,  the  integral  form 
in  A. 2  cannot  be  put  in  a  closed  form  by  this  method. 

B.   PARTIAL  DIFFERENTIAL  METHOD 

Return  to  the  integral  form  given  in  A. 2: 


2   Tr/2 

W(A)  =  (£)     /   9  cos  9  J  (A9)  df 

f  o 

o 


Define  new  equations:  CCA)  and  S(A) 


tt/2 

CU)    =      /      8    cos    9    J    (A9)    d9  (A. 6) 

o 
o 


tt/2 
S(A)    =      /      9    sin    9    J    (A9)    d9  (A. 7) 

o 


Differentiate   to    A   as    C'(A)    and   S'(A) 


,  tt/2      9 

C'(A)    =    ^r  CCA)    =         /      9      cos    9    J    '(A9)    d9  (A. 8) 

dA  o 

o 
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tt/2   ? 

S'(A)  =  ■£-  S(A)      /   9   sin  9  J  '(A9)  d9         (A. 9) 
dA  o 

o 


Multiply  both  A. 8  and  A. 9  with  A: 

tt/2       ? 

AC'(A)    =      /      9      cos    0    J    T(A9)    d(A9)  (A. 10) 

o 
o 


tt/2       9 

AS'(A)    =      /      9      sin    9    J    '(A9)    d(A9)  (A. 11) 

o 
o 


Using  Bessel  identity  J  ' ( A9 )d( A9 ) =dJ  (A9),  substitute 
&  J      o  o 

it    into   A. 10    and  A. 11: 


tt/2       9 

AC'(A)    =      /      8      cos    9    dJ    (A9)  (A. 12) 

o 
o 


tt/2       ? 

AS'(A)    =      /      0      cos    9    dJ    (A9)  (A. 13) 

o 
o 


Expand  both  integral  in  A. 12  and  A. 13  into  partial  integration 

?  tt/2  tt/2  9 

AC'(A)    =    9Z    cos    9    J    (A9)  -         /       (29cos9-9    sin9)J    (A9)df 

o  '  o 

o  o 


tt/2  ? 

=    -      /    (29    cos9-9      sin9)    J    (A9)    d9  (A. 14) 

o 
o 

9  tt/2  ? 

AS'(A)  =  9   sin9  J  (A9)      (29  sin9+9   cos9)J  (A9)d9 

o       '  o 

o 
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tt/2  ? 

=  I  J  (Ix)  -   /   (29  sinG+0^  cos6)J  (A0)d0     (A. 15) 
2   o  I  o 

o 


Rewrite  both  equation  in  A. 14  and  A. 15,  using  definition 
of  C  and  S.   It  will  end  up  with  the  following  pair  of  differ- 
ential equations: 

AC  (A)  =  -2CQ)  +  S'  (A)  (A. 16) 


AS'(A)  =  ~   J  (J\)  -  2S  -  C  (A. 17) 

l      o  l 


Rearrange  the  last  two  equations 


AC  +  2C  +  S'  =  0  (A. 18) 


XS'  +  2S  +  Cf  =  i  J  (5-X)  (A. 19) 

2   o  2 


The  last  two  pairs  of  equations  show  that  it  may  give  a  single 
equation  in  terms  of  S  or  C.   And  that  means  it  will  end  up  in 
integral  form  of  C  or  S.   The  conclusion  is,  the  integral 
form  in  A. 2  cannot  be  put  in  a  closed  form  by  this  method. 

C.   TABLE  OF  INTEGRATION 

Taylor  series  expands  any  analytic  function  f(x)  into  a 
power  series  of  x.  Consequently  cos  x  may  be  expanded  into 
a  power  series  of  x  : 
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cos  x  =  f(x2)  (A. 20) 


Substitute  A. 20  into  A. 2: 


f       2 
w(x)  =   /   x  f(x  )  J  (fx)dx  (A. 21) 

o  ° 


where  we  change  the  integration  variables  with  the  dummy 
variables  x  and  f. 

From  Table  of  Integration  [Ref.  10],  the  integral  form 
shown  below: 


v 
I0,(x)  =   /   tZK  L   J  (t)dt  (A. 22) 

ZK.  O 

O 


can  be  written  in  a  closed  form.   The  term  t(i  (t)  in  the 

Jo 

above  equation  is  replaced  using  Bessel  equation: 

tJ  (t)  =  -tJ  "(t)-J  '(t). 
o         o      o 

The  integral  form  in  A. 2  2  becomes: 


Iov(x)  =   /   t2k{-tJ  "(t)-J  '(t)}  dt 
2k  o      o 

o 


=  _  /  t2k+1J  '(t)dt-  /  t2kJ  '(t)dt 
o  o 

o  o 


=  -x2k+1J  »(x)+(2k+l)  /  t2kJ  '(t)dt-  /  t2kJ  '(t)dt 
o  o         _      o 

o  o 
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X 


Th 


=  -x2k+1J  '(x)+2k(x2kJ  (x)-2k   /  t2k-1J  (t)dt} 


I01  (x)  =  2kx2kJ  (x)-x2k+1J  '(x)-4k2I0,  -,(x)        (A. 23) 
2k  o  o  2k-l 


e  last  equation  shows  that  I„,  (x)  reduces  to  lower  order 


I  ,  , (x) .   This  form  is  called  reduction  form.   For  k=o, 
A. 2  2  and  A. 2  3  become: 


x 

I  (x)  =  /  t  J  (t)dt  =  =xJ  '(x)  (A. 24) 

o  o  o 

o 


The  equation  in  A. 24  agrees  with  the  property  of  Bessel 
function . 

If  the  reduction  form  applied  repeatedly  and  written  in 

a  general  form,  A. 2  3  may  be  written  as: 


x   ?1c 

f   t  *tJ  (t)dt  =  2k.Poly(2k).J  (x)+x.Poly ( 2k) . J  ' (x)    (A. 25) 

o  J  o  J  o 

o 


Lets    return   to    equation  A. 20 


00 

7  ?\c 

cos    x    =    f(x    )    =      E      fk   x  (A. 26) 

k=o 


Substitute  A. 26  into  A. 25: 


V  oo  oo 

/   f(t2)t   J    (t)dt    =      Z      a,  x2k    .    J    (x)+      E      b1x2k+1,(J    '  (x) 

o  ,  k  o  ,  k  o 

o  k=o  k=o 

(A. 27) 
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Apply  differentiation  to  both  left  and  right  hand  side  of  A. 27 


OO  00                                          00 

Z  f,x2k+1J  (x)  =  J  (x)  I  2ka,x2k~1  +  J  '(x)  Z   (2k+l)b,x2k 

,  k      o       o     i  _      k       o     ,           k 

k=o  k=o                 k=o 


+  J  '(x)  Z  a,x2k  +  J  "(x)  I    b,  x2k+1 
o        k       o        k 


Replace  J  "(x),  using  Bessel  equation 


J  "(x)  =  -  x_1J  '(k)  -  J  (x)  (A. 28) 


00  OO 

Z   f  x2k+1J  (x)  =  J  (x){  Z  2ka,  x2k-1  -  Z  b,x2k+1} 
k=Q   x      o       o      1     k  k 


+  J  »(x)  {Za,x2k  +  Z  (2k+l)b,  x2k-Zb,  x2k} 
ok  k      k 


00  00 

J  (x){  Z  (2k+l)a1 +.x2k+1-  Z  b1x2k+1} 
o  kt  k 

o  o 


00 

+  J  '  (x)  {E(a,  +2kb  )x2k} 
o        k    k 

o 


Equate  the  coeficient  of  the  equal  term  J   and  J  '  from  the 
left  hand  side  and  right  hand  side: 
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00 

00 

Z    fnx2k+1    =    Z    {(2k+l)aV4.n-bv}    x2k+1  (A. 29) 

ok  -,  k+1      k 


0  =  Z  (ak+2kbR)x  K  (A. 30) 


Therefore 


f   =  (2k+l)ak+1-bk  *  ak=   k-1^k-1  (A. 31) 


0  =  ak+2kbk       -  bk  =  -  ^  (A. 32) 


These  last  two  equations  give  the  solution  to  the  integral 
of  the  cos  x  form.   Lets  expand  cos  x  into  power  series: 


2       4       6 
cos  x  =  f   +  f,x   +  f„x   +  f 0x   + 
o     1       I  3 


then  using  relations  given  in  (A. 30)  and  (A. 31) 


a   =  0  b   =  0 

o  o 


f  f 

o  h   =    — 

'12  1    "  4 


fn-f  /4  f,-f  /4 

1   o  h        1   o 

a2  "     4  2  "  "    16 

f.-f  /4  f._-f  /4 

-   i   o  x   =     1   o 

3  "32  3        192 

etc . 
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The  conclusion  is,  we  succeeded  in  solving  Hankel  trans< 
form  of  cos  x  function,  but  we  failed  to  express  it  in  a 
closed  form. 
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